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Abstract 

Uniform-DFT filter banks are an important class of filter banks and their theory is well known. 
One notable characteristic is their very efficient implementation when using polyphase filters and the 
FFT. Separately, linear phase filter banks, i.e. filter banks in which the analysis filters have a linear 
phase are also an important class of filter banks and desired in many applications. Unfortunately, it 
has been proved that one cannot design critically-sampled, uniform-DFT, linear phase filter banks and 
achieve perfect reconstruction. In this paper, we present a least-squares solution to this problem and 
in addition prove that oversampled, uniform-DFT, linear phase filter banks (which are also useful in 
many applications) can be constructed for perfect reconstruction. Design examples are included which 
illustrate the methods. 
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I. Introduction 

In the last decade research in filter banks and subband processing has occurred at a fast pace [1*2]. 
[16], [17]. In particular, subband processing of signals is now common in a number of applications 
such as audio, image, and video encoders/decoders; acoustic echo cancelers used in teleconferencing 
systems; and spread spectrum communications systems [1]. In addition, emerging applications include 
multiple target tracking in radar, high-data-rate digital receivers for satellite channels, and suppression 
of interference signals in wireless applications [7], [13], [6], [2]. As compared to equivalent fullband 
processing, subband signal processing in these applications often leads to a performance increase due 
to the frequency band decomposition and potentially a reduction in computation due to the lower 
sampling rate of the subband (component) signals. 

An analysis bank decomposes (analyzes) a signal, x by partitioning its spectrum into frequency 

bands or subbands by using a set of M bandpass filters (BPFs) denoted h 0 h\/_! in Fig. 1. Once 

partitioned, the signals are then downsampled by a factor of D (due to their reduced bandwidth). 
Upon analysis, each of these subband signals (tq, may be processed (subband processing). 

After processing, the subband signals are synthesized into the fullband signal. ,r, by first upsampling 
(zero insertion) to the original sampling rate followed by bandpass filtering (fo, .... f\/_i) to remove 
spectral images introduced from upsampling. Such a system of upsamplers and bandpass (synthesis) 
filters is referred to as a synthesis bank. The combination of analysis and synthesis banks is called a 
filter bank [17]. 
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II. Untform-DFT Filter Banks 

One of the more popular designs for filter banks, called the uniform-DFT filter bank, assumes that 
all analysis, synthesis filters are frequency-shifted versions of an analysis, synthesis prototype lowpass 
filter (LPF), ho, fo respectively [17]. The time-domain relations are then given by 

h m (n) = h Q {n)e^ m ^- (N - x]l2]IM 

fm(n) = (1) 

where m = 1, . . ., M - 1 denotes the subband and N is the length of the filter. The frequency-domain 
relations are then given by 

H m (e JW ) = Ho (e j(w_2ffm / A/) ) e -jrrn(JV ' _1) / A/ 

F m (e j “) = F 0 (e j( ^~ 2ffm/A7) ) (2) 

A typical magnitude response of analysis/synthesis filters related by (2) with M = 4 is given in 
Fig. 2. 

A. Polyphase , Uniform-DFT Filter Banks 

In the analysis bank shown in Fig. 1, much computation is wasted in the implementation since 
after filtering, only one out of every D samples is retained. Similar filtering inefficiencies occur in 
the synthesis bank since only one out of every D samples is non-zero prior to synthesis filtering. A 
more efficient but equivalent implementation, called a polyphase, uniform-DFT filter bank relies on a 
polyphase representation of h 0 and f 0 (prototype LPFs) and a discrete Fourier transform (DFT) and 
inverse DFT (IDFT) as shown in Fig. 3 [12], [17]. Note that in Fig. 3, / — M/ D is referred to as the 
oversampling factor and will be discussed in the next section. The polyphase filter coefficients used in 
Fig. 3 are given bv 

e m (w) = ho(nD - m) 

r m {n) = f 0 {nD - m) (3) 

where m = 0, . . . , M - 1. As an example, if we assume M = 4, D — 2. and a causal. FIR analysis filter 
prototype, ho the polyphase filters would be given by 

eo = [M0),M2),M4)./>o(G),..f 
ei = [0,Ml),M3),M5),...f 
e 2 = [0,M0).M2),/>o(4),...f 

e 3 = [0, 0, fio(l), /?o(3 ), . . ] T . (-1) 
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With the polyphase representation, analysis and synthesis filtering occurs at the lower, subband sam- 
pling rate resulting in a lower processing speed. In addition, the DFT and IDFT are implemented 
with a fast-Fourier transform (FFT). These combined effects result in a significant computational 
reduction. 

In most of the literature, filter bank design is aimed at perfectly reconstructing x (to within a scale 
factor and delay) from the critically-sampled ( D = M or I = 1) subband signals x*o, (Fig. 3). 

Critically-sampled refers to the fact that the number of fullband samples per second equals the total 
number of subband samples per second. In this case, one must carefully control the aliasing in the 
subband signals through proper design of the analysis and synthesis filters. Many critically-sampled 
designs are available that yield the perfect reconstruction property [15]. 

B. Oversampled , Polyphase , Uniform- DFT Filter Banks 

A fundamental problem in processing critically-sampled subband signals is the large level of subband 
aliasing which is present regardless of the analysis and synthesis filter designs due to the overlapping 
filters. With such a large amount of aliasing (which acts as a disruptive noise signal), many subband 
signal processing applications do not perform well and may require a large additional overhead to 
partially compensate for the aliasing [5]. 

One solution to the subband aliasing problem in applications is to oversample ( D < M) the subband 
signals to minimize the effect of the aliasing [12], [3]. We typically choose the oversampling factor to be 
close to unity in order to minimize the data rate or choose 1 — 2 for simple implementation as in Fig. 3. 
In the latter case, this implies twice the total number of samples per second in the subbands than 
are required at the fullband rate. Howver, even with oversampled subbands, many signal processing 
applications can have fewer total computations per sample as compared to the fullband case [4]. 

III. Linear Phase Filter Banks 

Linear phase filter banks are designed so that the analysis filters have a linear phase response. 
This should not be confused with the idea that the frequency response of the filter bank has linear 
phase, i.e. that the end-to-end impulse response of the filter bank is symmetric. There are many filter 
bank applications that require linear phase analysis filters such as subband image coders and filter 
bank-based digital receivers [13], [6], [18]. The basic theory of linear phase filter banks with perfect 
reconstruction is also well-known (for complete details, the reader is referred to [11]). 
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IV. Uniform-DFT, Linear Phase Filter Banks 

It has been shown that critically-sampled uniform-DFT filter banks in which the analysis filters 
have linear phase cannot be designed to have perfect reconstruction [8], [9]. In this section, we analyze 
the uniform-DFT, linear phase filter bank and formulate a least-squares synthesis filter solution in 
the critically-sampled case which approximates perfect reconstruction. In the oversampled case, we 
formulate a synthesis filter design that leads to perfect reconstruction. 

.4. Time-Domain Analysis 

In order to formulate the design for the uniform-DFT linear phase filter bank, we begin by charac- 
terizing the family of responses to a unit-pulse with delay. I or S(n — l) for 1 </<£). In general, the 
filter bank is not time-invariant but rather periodically time-invariant with period D thus resulting in 
our need to consider this family. The subband signals (Fig. 1) in response to S(n - l) can be written 
as 

After upsampling and synthesis filtering the subband signals in (5), we have 

t/m,/(n) = 52h m (kD-l)f m (n-kD) (6) 

k 

as in Fig. 1. By substituting the prototype filter relations in (1) into (6) we have 

y mJ (n) = £ ho(kD - l)fo(n - (7) 

k 

The response of the filter bank to 8{n - l) is then sum of the signals in (7) or 
A/-1 

i ‘( n ) = T y™A n ) 

m= 0 

l _ gj2jr[n— /— ( JV— 1 )] 

= ho(kD - l)f 0 (n - kD) _ j 2 7r [n-/-{.v- 1 )]/A/ 
k 

= M L h ^ kD - 0 /o l l + (* v “ !) - kD - TV] & [n - / - (A r - 1) + <jM] (8) 

9 k 

since 

l _ € J‘2r[n-l-{.\-l)] 

r ‘V-D+.A/J. <»> 

c 9 

One condition for perfect reconstruction (with delay) is that the family of responses to the delayed 
unit-pulse is equal to a family of equivalently delayed unit-pulses. Specifically we require 

xi(n) — S[n-l-(X- 1)] (10) 
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for 1 < l < D. 

We can arrange the conditions in (10) using (8) in matrix form as 


Hfn = 


where 


h Q {D + 1) 


ho(D) 

0 


0 

*o(l) 




*o(0) 

0 


0 ho(<?max, 1 £) 0 *o(<? max,] D-D) 

*o(W*, 2 £>+ D - L) 0 MWx.20-0 0 


with dimension 


dim(H) 


r D 

L/=i 


N - 1 + l 


M 


- 1 


(11) 


0 


RT 


o o hoiqmax.D- 1 ^ ' 

>»o(W.x.£>£ + D 0 0 


X A r , 


(13) 


•jfmax.l — 


;v - 1 + i 

D 


(14) 


and 


£?. Critically-Sampled Case 


E,Tl i T7 ±i ]- D ENl^J- 1 


i T 


(15) 


In the critically-sampled case, D = M we consider the matrix form of (10) for each / individually. 
This is equivalent to decomposing the matrix equation in (11) into a series of matrix equations [due 
to the decoupled nature of (11)] as 


for 1 < l < D where 


h 0 (M - l) 
h 0 (2M - l) 


H;fj = yS, 


0 

h 0 (M - l ) 


0 

0 


(16) 


0 

0 


H, 


*0( Qmax.l /) /lo [(^fmax.i l)il/ /] 


h 0 (2M - !) 


h 0 (M-l) 


( 


0 

0 


o 0 hoUjmux.iM — l) /io[(w - 1)4/ - /] 

0 0 0 h 0 (q max ,iM - I ) 
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f i 


fo [(^ — 1) + / - ?ma xjM] 
fo [(N - 1) + / - (W,/ - 1) A/] 

/ 0 [(;V-1) + /-2A/] 

/o[(iV-l) + /-A/] 


(IS) 


and 

Note that only a subset of the synthesis prototype coefficients are present in (18). However, over the 
full range of l all coefficients are accounted for. 

By choosing the filter length, N = pM for integer p, the matrix H/ is (2 p - 1) x p and f/ is p x 1. 
For this case, our problem is overdetermined and as well known, there is no solution for the set of 
fy , f 2 • • ■ - f\/ that constitute the synthesis prototype filter. f 0 and leads to perfect reconstruction [9]. 
However, we can formulate a new, least-squares (LS) design problem as follows. 


0 1 0 


(19) 


Given: H/ € ft (2p_1)xp , Si G £ (2p ~ 1)xl for 1 < / < M 
Minimize: V'(f/) = ||<5; — H;f ;|| 2 for 1 < / < M. 


The LS solution to this problem is 


f( = H + 8 , 


( 20 ) 


where 

H+ = (H”H/) _1 (21) 

is the pseudoinverse of H/, det(H*H/) ^ 0, and H“ is the complex conjugate transpose of H/. 

For the choice N = pM + q where 1 < q < M - 1, f/ will have a different length according to I 
but the problem remains overdetermined and the least-squares solution in (20) still holds. General LS 
approximations of perfect reconstruction filter banks are discussed in [14]. 

By design, the analysis prototype LPF, h 0 has a linear phase and symmetric impulse response of 
the form 


h 0 (n) = h 0 [(A r -l)-n] 


where n = 0 A’ — 1 . From (1) it is clear that the analysis filters, hj , . . . , h v-i will also have linear 

phase responses. We can prove that under the formulation in (11), the LS solution for the synthesis 
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prototype LPF fo will also have a linear phase and hence the synthesis filters will also 

have linear phase responses. The proof is as follows. 


Consider the case where N = pM. We have q m&x - p from (14). and applying (22) to (17) and 
substituting the result into (16) we have 


h 0 [(N -\) + l- M] 0 0 


" 0 " 

h 0 [(N - 1) + l - 2M] : 


fo(l-l) 



: 0 


M-l + M) 


0 

h 0 (l-l) h 0 [(N -l) + l-M] 



= 

1 

0 h 0 [(N -\) + l-2M] 


f 0 [(N- 1) + /-2.U] 


0 



. foUN-l) + l-M] . 



0 0 ho(l-l) 


_ 0 _ 


Now consider the vector 


fo (M - 1) 
fo (2 M - /) 

fo (-V - M - /) 

fo(N-l) 


(24) 


and the equivalent to (23) which uses f) instead of f) 


1 

o 

1 

o 

0 


0 

h 0 (N-M-l) ••• 



fo(U-i) 




0 


fo (24/ - /) 


0 

h 0 (M - l ) 

h 0 (N - l) 



= 

1 

0 

h 0 (N-M-l) 


fo (A 7 - M - I) 


0 

• 



fo (A’ - /) 



0 

0 h 0 (M - l) 


0 


Equation (25) is identical to (16) except that elements of the matrix along any diagonal are in reverse 
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order. We can rewrite (25) as 


1 

O 

£; 

i 

o 

0 


0 

h 0 (‘2M - l) ••• 



fo(N ~ l) 


\ 

* • , 

0 


/o(-V - M - l) 


0 

ho(N - l) 

ho(M-l) 


• 

~ 

1 

0 

••• h 0 (-2M -l) 


fo (2AJ - /) 


0 




fo(M-l) 



0 

0 ho(N-l) 


0 


Using the LS solution given in (20) in both (23) and (26) yields fo that satisfies 

/o(«) = fo(N-l-n) 


or in other words, fo also has linear phase. The same result can also be shown for the case where 


N = pM + q where q = 1 , . . . , M — 1 . 

C. Oversampled Case with Integer Oversampling Factor 

We next consider the oversampled case where the oversampling factor. / > 2 is an integer. In this 
case, we use (16) but with 

ho(D-l) 0 0 

h 0 (2D-l) ho(D-l) 0 0 


Hi 


ho(<lmaxJ D l) 


ho [Omax.l - 1 )D~l] 


ho(‘2D-l) 


ho(D-l) 


(28) 


0 

0 


0 


ho{q max, iD-l) 
0 


ho [(^maxJ ” \)D- l] 

ho (^niaxj I-) ~~ l) 


fo [(A” ” 1) + / — 7m ax,/ D\ 
fo [(jV — 1) -f / — (7max./ ” 1) D\ 


1 1 


(29) 


fo [(A ,r — 1) -I- / — 2D] 
fo[(N- 1)+Z-D] 


l = 1, . . . , D, and is given by (14). 

By choosing N = pM for integer p, the matrix H/ is (2 p - 1) x Ip and f ] is Ip X 1. For this case, 
our problem is underdetermined and there exists an infinite number of solutions for f/ and hence fo 


December 2. 1999 


DRAFT 



10 


SUBMITTED OCT. 13. 1999 TO IEEE TRANS. SP 


that lead to perfect reconstruction. We therefore reformulate the design problem so as to choose the 
fo that has minimum stopband energy. Assume the cutoff frequency of the analysis prototype filter is 
u c — 7 t/A/. Then the synthesis prototype LPF we choose minimizes the cost function 


where 


J(io) = f \F 0 {e^)\ 2 d^ 

JlMr 


jV-1 

F 0 (e j n = ^/o(n)e"^". 

n=0 


Substituting (31) into (30) we have 

A r -1 r 


-'(fo) = Y 


«(») 


n = 0 L 


M - 1) 

~i : 


2tt 


A'-l 


m = 0 


A’-l- 

E 

i 


E /o(m) E /o(m + fc)sinc 


(A/ - 1)jt t 2;r t 

to - 77*0 sto 


A/ 


A/ 


where sine (a:) = S1 ^ X ' and 


isinc(^) 


^sinc(^) ±sinc(-^) 
0 |sinc(j^) 


isineC^) 

2 s * nc ( ^\7^) 


fSincf^p) isinc(Nr) 


^sinc(^j) 0 


We rewrite (32) as 


./(fo) = ELJ> f r ufo 


(30) 


(31) 


(32) 


(33) 


(34) 


where 


U = I - 


-)S 


(35) 


M - 1 

is a real, symmetric matrix. Since J(fo) > 0 for all fo, U is a positive definite matrix and there exist 
a non-singular matrix C that satisfies 


c r uc 


Ao 


Ajv-i 


(36) 


where A 0 , • - •, A.v_i are the positive eigenvalues of U. Define 


fo = Pf 0 


(37) 
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where 


\Ao 


P = 


i-i 




, r -i 


Using (36)-(38) in (34) we have 


T(f) _ 

J (to) - * 0 * 0 * 


Substituting f 0 = P L fo into (11) vve have 


1 


HP- A fo = —A 


or 


Hf „ = jja 


where 


H = HP' 


(38) 


(39) 


(40) 


(41) 


(42) 


and has dimension — D) x N . We can now formulate the design problem for the oversampled case 
(integer /) under the minimum stopband constraint. 

Given: H 6 ft (2 r -D)X;V , A 6 

Minimize: V’(fo) = ||fo|| 2 
Subject to: Hfo = -^yA. 

The LS solution to this problem is 

fo = (43) 

where 

H + = H' (Mj" 1 (44) 

is the pseudoinverse of H, det(HH~) ^ 0. and H~ is the complex conjugate transpose of H. Given 
the solution in (43), the synthesis prototype filter coefficients are then given by 

fo = P _1 fo- (40) 
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D. Oversampled Case with Rational Oversampling Factor 

The dimension of H given in (13) which is a function of AT, A/, and D determines whether the design 
problem is over-determined (no solution for Fq) or under-determined (an infinite number of solutions 
for fo). Depending on the case, we may design f 0 according to the previous cases. 

V. Design Examples and Performance Analysis 
The design algorithm for uniform-DFT linear phase filter banks is as follows. 

1. Design an analysis prototype LPF, ho with linear phase 

2. (Critically sampled case) Determine H / for each / as in (17) 

(Oversampled case) Determine H as in (42) 

3. (Critically sampled case) Compute f/ for each l according to (20) 

(Oversampled case) Compute fo according to (45) 

4. The synthesis prototype LPF, fo is constructed according to (18) or (45) for the critically sampled 
or oversampled case, respectively, fo will be linear phase. 

A MATLAB program which implements the above algorithm is available at [10]. We present several 
design examples illustrating the method. 

Example 1 (Critically-sampled) Here we have M = D — 4 and the analysis prototype LPF is de- 
signed with a length N = 32, Hamming-windowed sine function (firl in MATLAB). Figs. 4 and 5 
illustrate magnitude responses of the prototype filters. Fig. 6 illustrates the worst-case impulse re- 
sponse, i.e. the response to S(n — l) that has the largest magnitude component (artifact) other than the 
main, delayed impulse. Fig. 7 illustrates the magnitude response of the filter bank for the worst-case 
impulse response. 

Example 2 (Oversampled) Here we have M — 4. D — 2 and the analysis prototype LPF is designed 
with a length N = 32, Hamming-windowed sine function (firl in MATLAB). Figs. 4 and 8 illustrate 
the magnitude responses of prototype filters. Fig. 9 illustrates the impulse response, and Fig. 10 
illustrates the magnitude response of the filter bank for the worst-case impulse response. 

VI. Conclusions 

This paper has presented a least-squares solution to the design of uniform-DFT linear phase filter 
banks. In the critically sampled case, the design yields near perfect reconstruction while in the 
oversampled case, the design yields perfect reconstruction. An algorithm and examples were presented 
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to demonstrate the method. 


Ac K N O VV L E D G M ENTS 

The authors would like to acknowledge the useful discussions with Dr. Dennis Morgan of Bell 
Laboratories, Lucent Technologies and Dr. Henrique Malvar of Microsoft C’orp. 

References 

[1] A. A kans u and M. Smith, Subband and Wavelet Transforms: Design and Applications , Norwell. MA: Kluwer Academic 
Publishers, 1996. 

[2] A. Akansu, M. Tazebay, M. Medley, and P. Das, “Wavelet and subband transforms: fundamentals and communication 
applications/ 1 IEEE Comm . Mag., pp. 104-115, Dec. 1997. 

[3] Z. Cvetkovic and M. Vetterli, “Oversampled filter banks, 11 IEEE Trans. Signal Processing . vol. 46. pp. 1245-1255. 
May 1998. 

[4] P. De Leon, “Optimization of the LMS subband, adaptive filter system, 11 in Proc. 1997 International Conference on 
Signal Processing Applications & Technologies , 1997. 

[5] A. Gilloire and M. Vetterli, “Adaptive filtering in subbands with critical sampling: analysis, experiments, and 
applications to acoustic echo cancellation, 11 IEEE Trans . Signal Processing, vol. 40, pp. 1862-1875, Aug. 1992. 

[6] A. Gray, P. Ghuman, and S. Sheikh, “High rate digital demodulator ASIC, 11 in Proc. 1998 International Conference 
on Signal Processing Applications & Technologies . 1998. 

[7] L. Hong, “Multiresolutional multiple-model target tracking, 11 IEEE Trans, on At rasp. Electron. Syst vol. 30. pp. 518- 
524, Apr. 1994. 

[8] H. Malvar, “Modulated QMF filter banks with perfect reconstruction. 11 Electron. Lett., vol. 26. no. 13. pp. 906-907. 
Jun. 21, 1990. 

[9] H. Malvar, Signal Processing with Lapped Transforms , Norwood, MA: Artech House, 1992. 

[10] http : //telsat .mnsu.edu/~pdeleon/uniDFT_linph_fb.rn 

[11] T. Nguyen and P.P. Vaidyanathan, “Structures for M-channel perfect reconstruction FIR QMF banks which yield 
linear-phase analysis filters, 11 IEEE Trans. Acoust. Speech Signal Processing, vol. 38, no. 3, pp. 433-445. Mar. 1990. 

[12] L. Rabiner and R. Crochier, Multirate Digital Signal Processing , Englewood Cliffs, NJ: Prentice-Hall. 1983. 

[13] R. Sadr, P.P. Vaidyanathan, and D. Raphaeli, and S. Hinedi, “Parallel digital modem using multirate digital filter 
banks, 11 JPL Publication , 94-20, Aug. 1994. 

[14] M. Spurbeck and C. Mullis, “Least squares approximations of perfect reconstruction filter banks." IEEE Trans. 
Signal Processing, vol. 46, pp. 968-978, Apr. 1998. 

[15] K. Swaminathan and P.P. Vaidyanathan, “Theory and design of uniform DFT, parallel, quadrature mirror filter 
banks, 11 IEEE Trans. Circuits Syst., vol. CAS-33. no. 12, pp. 1170-1190. Dec. 1986. 

[16] P. P. Vaidyanathan, * ’Multirate digital filters, filter banks, polyphase networks, and applications: a tutorial." Proc. 
IEEE , vol. 78, no. 1, pp. 56-93, Jan. 1990. 

[17] P. P. Vaidyanathan, Multirate Systems and Filter Banks, Englewood Cliffs, NJ: Prentice- Hall. 1993. 

[18] J. Woods and S. O’Neil, “Subband coding of images,” IEEE Trans. Acoust. Speech. Signal Processing, vol. ASSP-34. 
pp. 1278-1288, Oct. 1986. 


December 2, 1999 


DRAFT 















Magnitude (dB) Magnitude (dB) 


SUN AND DE LEON: DESIGN TECHNIQUES FOR 



December 2, 1999 









December 2 . 1999 




Magnitude (dB) 



DRAFT 


December 2. 1099 






